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ABSTRACT 

Accretion discs are common in binary systems, and they are often found to be misaligned 
with respect to the binary orbit. The gravitational torque from a companion induces nodal 
precession in misaligned disc orbits. We calculate whether this precession is strong enough to 
overcome the internal disc torques communicating angular momentum. For typical parameters 
precession wins: the disc breaks into distinct planes that precess effectively independently. 
We run hydrodynamical simulations to check these results, and confirm that disc breaking is 
widespread and generally enhances accretion on to the central object. This applies in many 
cases of astrophysical accretion, e.g. supermassive black hole binaries and X-ray binaries. 
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1 INTRODUCTION 

Accretion disc s (e.g. IPrinde & Rees 197^ : IPrinde 198II : 
iFrank et al. 20021) appear in many astrophysieal systems. In 
most eases these discs are probably not eompletely axisymmetrie. 
Discs may be locally tilt ed by the Lense-Thirring eff ect of a central 
misaligned blaek hole (iBardeen & Petterson. 1975h . by radiation 
dPringle 199dPringle 19971) or by the gravity of a eompanion (e.g. 
iLubow & Qgilvie. 20001) . Discs in or around supermassive black 
hole (SMBH) binaries formed by galaxy mergers may be mis¬ 
aligned with respect t o the binary orbit through the chaotic natur e 
of AGN aceretion (iKing & Pringle 200dKing & Pringle 200^ . 
The effeet on the dise in all these cases is similar. The lack of 
symmetry produees a torque on misaligned rings of gas whieh 
makes their orbits precess differentially. Given a sufficiently 
strong viscosity communicating the precessio n between the rings, 
the disc warps. IPanaloizou & Pringle (1983h showed that warps 
ean propagate in two distinct regimes: wave-like for a < H/R, 
and diffusive for a > H/R where a is the Shakura & Sunyaev 
dimensionless viscosity parameter (IShakura & Sunvaev. 1973h 
and H/R is the disc angular semi-thickness. In this paper we foeus 
on diffusive systems with a > H/R, which typically holds for 
accretion discs around black holes. 

For diffusive discs subject to differential preeession, the ex¬ 
pected evolution is that dissipation through viscosity allows the 
inner parts of the disc to align, joined by a smoothly warped re¬ 
gion to the still misaligned outer parts. In the case of a misaligned 
disc around a spinning black hole, this is often called the Bardeen- 


* suzan.dogan@ege.edu.tr 
t chris.nixon@jila.colorado.edu 
$ Einstein Fellow 


Petterson effect. Until recently it was implicitly assumed that this 
is what always occurs in a diffusive disc, i.e. that the internal disc 
torques would always be able to eommun icate the preeession. How¬ 
ever an analytic study bv lOgilvie (1999h pointed out that the effect 
of a disc warp was to weaken the communication of angular mo¬ 
mentum in the dise, and so weaken the disc’s ability to hold itself 
together. Although his study assumed a locally isotropic viscous 
process, there appears no reason to assume this behaviour does 
not hold for a viseosity drive n by turbulenee, sueh as the magne- 
torotational instability fMRI: iBalbus & Hawley. 199lh . Indeed for 
a viscosity driven by magnetic fields it is likely that the vertical 
viscosity, associated with keeping the disc flat, is weaker than that 
of an isotropic model, as the vertical gas shear is probably oscil¬ 
latory whereas the azi muthal shear gr ows secularly as gas parcels 
continually move apart (IPringle. 1992h . 


We note that it is conventional to use the term isotropic vis¬ 
cosity, but that this can be misleading. For a warped disc this term 
means that the horizontal and vertical shear in the warp are as¬ 
sumed t o be damped by viscous dissipation at the same average 
ra tes (cf. Lodato & Pringle. 20071 . e g. 40). This as sumption, made 
bv IPapaloizou & Pringle (1983h and lOgilvie fl999h . leads to the re¬ 
sult that the azimuthal shear viscosity vi oc a, but that the vertical 
viscosity V 2 oc I/a, quite eontrary to any naive belief that vi and 
V 2 might end up roughly equal. This apparently paradoxical result 
eomes about beeause a small value of isotropie viscosity a allows 
a large resonant radial velocity vr in the warp: the viscous dissi¬ 
pation rate goes as av^ oc 1/a and so increases as a deereases. 
A warped-disc isotropic viscosity explicitly does not assume that 
V] = y 2 , as was the case in the early work on warped discs by e.g. 
iBardeen & Petterson fl975h . 


An a viscosity that acts isotropically, as described above, ap¬ 
pears to hold for viscosity arising from turbulence induced by the 
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Figure 1. 3D surface rendering of the warped disc after 0 and 6.5 binary 
orbits. The disc was initially inclined at 10° to the binary plane with no 
warp. The disc is viewed along the binary orbital plane and the arrow points 
in the direction of the companion. 


MRI (iKing et al.. 2013h . However, very little effort has been di¬ 
rected towar ds estimating the effecti ve viscosities in magnetised 
warped discs. iTorkelsson et al. (200Qh performed shearing box cal¬ 
culations to follow the decay of an imposed epicyclic shearing 
motion, which mimics a warp. Their results are in approximate 
agreement with an isotropi c viscosit y (which predi cts a 2 ~ 1/cr; 
IPapaloizou & Pringle 1983h . Further. lOgil vie f200l) developed an 
analytic model for the dynamical evolution of magnetorotational 
tur bulent stresses which p redicts agreement with the conclusions 
of ITorkelsson et al. (20001) . Both of these investigations, numer¬ 
ical and analytic, allow the effective viscosity from magneto¬ 
hydrodynamic (MHD) turbulence to be anisotropic, but both con¬ 
clude that it is close to isotropic. 

The realization that the viscosity may not be strong enough to 
hold the disc together has significant implications. If the viscosity 
is too weak, or the external torque on the disc too strong, the disc 
may instead b reak into distinct plane s with only tenuous gas flows 
between them (iNixon & King. 20 13) . If in addition these planes are 
sufficiently inclined to the axis of precession, they can precess until 
they are partially counterrotating, promoting angular momentum 
cancellation and rapid infall - disc tearing. Tearing occurs in discs 
inclined to the spin of a central black hole (iNixon et al.. 2012h and 
in a circumbinary d isc around a misaligned central binary system 
jNixon et al.. 2013h . 

In this paper we want to And out if tearing can happen in¬ 
side a binary, i.e. if a disc around one component can be disrupted 
by the perturbation from a companion. This would have signifi¬ 
cant implications for all binary systems: e.g. f uelling SMBH dur - 
ing the SMBH binary phase (cf. Figs. 6 & 7 of INixon et al. 2013b 
and accretion outbursts in X-ray binaries. In the Lense-Thirring and 
circumbinary disc cases, disc breaking starts from the inside and 
works its way outwards. But we also want to know if breaking and 
tearing can instead start from the outer edge of a disc internal to 
a binary, and work its way inwards. To do this we consider binary 
systems with an initially planar disc around one component, mis¬ 
aligned with respect to the (circular) binary orbit. Following the 
methods of our earlier papers on disc tearing, we first compare the 
disc precession torque with the disc viscous torque to determine 
whether the disc should warp or break. Then we check our findings 
by comparing this result with hydrodynamical simulations. 


2 TEARING UP THE DISC 


The disc precession caused by th e presence of a bi nary companion 
is retrograde, and has frequency (iBate et al.. 20001) 

ilp = (-] ilcose. (1) 

Here 6 is the inclination angle between the disc and the binary, Mq 
& Ml are the masses of each component of the binary with the disc 
around Mi, a is the binary separation, R (assumed a) is the disc 
radius, and Q - (GMilR^y^^ is the disc orbital frequency. 

We get an idea of whether the disc tears by estimating the 
disc precession frequency in a typical case. The disc cannot extend 
past the Roche lobe r adius (more preci sely the tidal truncation ra¬ 
dius Rtide ~ 0.87 Rpt : [Frank et al. 2002h . so we take Mi = Mq and 
^max ~ 0.35a. Putting this into Q gives 


^^p,max ~ 0.03fl cos 6. 


( 2 ) 


So in this case the precession time is only ~ 30 dynamical times, 
suggesting that tearing is possible, as the viscous communication 
in the disc is likely to be significantly slower than this. 

We expect the disc to break when the precession induced in 
the disc is stronger than any internal communication in the disc. 
This communication can be due to the usual planar disc viscosity 
(vi), the viscosity arising from vertical shear in a warped disc (y 2 ) 
or pressure. In the simulations presented here we are focusing on 
the regime with a > HjR and therefore the communication due to 
pressure is small and we return to this point in Section|4] 

W e can write the magnitudes o f the viscous torques per unit 
area as jPapaloizou & Pringle. 1983h 


I I _ ?>nv{LR^Q. 
I ” InRH 


( 3 ) 


and 
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Here E is the disc surface density and / is the unit 
angular momentum vec tor. For Keplerian rotation and a 
IShakura & Sunvaev I1973h viscosity v/ = we can write the 

total as 


|Gto>all = |Gv,| + |G 


V2 


H 

2 'r 


[Sori + a2 |iAI] 


(5) 


where li/'l is th e warp amplitude and defined as \>jj\ = R\dlldR\ 
dOgilvie. 1999h . 

We can compare this to the magnitude of the precession torque 
per unit area 


|Gp 


l^p X L\ 


AMi\a) 


1,R^Q:^ cost? sin 0 


( 6 ) 


to give an idea of where in the disc we expect breaking to occur. 
Here L is the angular momentum density vector. To break the disc 
the precession must be stronger than its viscous communication, 
i.e. |Gp| ^ IGtotail, giving 


-^break ~ 


4 (ofi + f fj 

sin 2^ R Mo 


1/3 


a. 


(7) 


This break radius accounts for both the azimuthal and vertical vis- 
cositi es in a warped disc. In contrast, the p revious disc tearing pa¬ 
pers (INixon et al.. 2012lNixon et al.. 20131) used ai = a and con¬ 
sidered the initial conditions of a flat disc with \i^\ = 0. 
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Figure 2. 3D surface rendering of the disc which was initially inclined at 30° to the binary plane with no warp. These snapshots are taken after 0, 6.5, 17.5 and 
25.5 binary orbits. The disc is viewed along the binary orbital plane and the arrow points in the direction of the companion. In this simulation the disc breaks 
into two distinct planes after ~8 binary orbits. 


It is not straightforward to evaluate 0 as both a I 
and Qf 2 are strong func tions of the warp amplitude \if/\ 
dOgilvie. 199SlOgilvie. 200C)h and the warp amplitude itself is un¬ 
known before performing a full 3D calculation of the disc evolu¬ 
tion. In previous work it has sufficed to conservatively use ai = a, 
but to exclude the a 2 term. For large cr > 0.1 this is reasonable, but 
for smaller a the vertical viscosity is clearly important. Proceeding 
with the method of the earlier papers we get 


/ 4a H MA 
(sin26» « 


( 8 ) 


but we caution that this equation is not relevant for cr 0.1 and 
small inclination angles where the strong vertical viscosity can re¬ 
sult in rapid disc alignment. In such cases the total internal torque 
must be considered 0, but we note that this is not trivial to evaluate 
beforehand. 

We can evaluate 0 for typical parameters, giving 


( a \ 

^break ^ 0.16a — j 


a IhIrV'^ 


0.01 


Mo 


(sin20r‘/V 


(9) 


This disc tearing criterion is equivalent to requiring a minimum 
inclination of the disc to the binary orbit, ^min, defined by 


H Ml 

y- 

R Mo 



( 10 ) 


^break 


> 2.5 — 


Ml 


1/3 


where M - M\ + Mq is the total binary mass, so 
(H/R\ 


sin 26 > 0.06 


/ Of \ / H/R 

lo.iMo.oi 


(Ml 


( 11 ) 


cni) becomes 

< Mo) (12) 


since M ^ Mq in this case. 

If instead the disc is around the more massive binary compo- 
nent we have Mi > Mp and the disc size is approximately 0.6a 
(lArtvmowicz & Lubow. 1994i) . In this case, breaking occurs if 




HIR\lMilMo 
10 


(Ml > Mo) 


(13) 


For typical black hole disc parameters a - Q.l, HjR < 10"^ 
almost all discs with a reasonable misalignment should break (cf. 
[T2l[T3l). However, a very large mass ratio Mi/Mo » 1 makes the 
perturbation by the smaller black hole so weak that breaking would 
occur only after a very long interval. 

For X-ray binaries breaking can clearly be avoided in some 
cases, but probably occurs in others. First, if mass is trans¬ 
ferred by Roche lobe overfiow, the accretion disc forms closely 
aligned to the binary plane. So to get any disc inclination to 
the binary plan^i] in a close stellar-mass binary requires a torque 
to tilt the disc out of the plane. Here Her X-1 is interest¬ 
ing, as this system is known to have a tilted precessing disc 


We can simplify this formula in two limits. If the disc is around the 
less massive component we have Mi < Mq and the tidal limit on 
the disc size requires 


^ Note that the binary orbital plane and the spin plane of the black hole may 
be misaligned and so a disc a ligned to the binary plane could still experience 
Lense-Thirring tearing (e.g. lNixon & Salvesen. 20l4) . 
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Figure 3. 3D surface rendering of the disc which was initially inclined at 45° to the binary plane with no warp. These snapshots are taken after 0, 6, 20.5 and 
31 binary orbits. The disc is viewed along the binary orbital plane and the arrow points in the direction of the companion. In this simulation the disc initially 
breaks into two distinct planes after ~1 binary orbits. Then a third ring is broken off, but quickly interacts with the outer ring. 


(iKatz. 1973h. which sets limits on the viscosity coefficient a 
(cf. iKing et al.. 20131). The disc t i lt is plausibly attributed to ra¬ 
diation warping (Petterson 1977bl: IPetterson 19773 : IPringle 19^ : 
IWiiers & Pringle 19991 : lOgilvie & Dubus 200lh ^vided tha t the 
mass input occurs at small disc radii. Iwiiers & Pringle fl999h esti¬ 
mate a ^ 03,HIR ^ 0.04, and Mi/Mq ^ 0.5,R/a ^ 0.24. From 
cni) these give the requirement sin2^ > 1.7 for breaking to occur, 
which is of course impossible. This is reassuring, as the disc in Her 
X-1 appears to process as a single body. However, slightly larger or 
thinner discs, or ones with lower viscosity, could easily have values 
of sin 26 allowing for disc breaking. 


We note that /^break is the radius outside which we expect 
the disc to break. This is the opposi te of the Lense-T hirring 
jNixon et al.. 2012h and circumbinary (iNixon et al.. 2013h cases, 
and raises new possibilities. If for example a disc ring broken from 
the outer edge contained more angular momentum than everything 
inside, the outer disc might be able to sweep the entire inner disc 
in to the accretor and leave behind a single misaligned ring. This 


would presumably spread viscously and possibly repeat the pro¬ 
cess. 


3 SIMULATIONS 


To check our analytical reasoning above, we perform 3D hy- 
drodynamical numerical simul ations using the phantom smoothed 
particle hydrodynamics c ode (Mce. 2012h. as in previous stud- 
ies of warped discs (e.g. iLodato & Price 20 1Q|: iNixon 2012h and 
broken discs jNixon et al. 20121 : iNixon et al. 2013li . Disc break- 
ing has also been observed in the circumbinary simulations of 
iLarwood & Panalo izou (199 %. an d the unforced warped-disc sim¬ 
ulations of iLodato & Price f201Clh. _ 

We follow the method of Nixon et al. (2013h . but simulate 
discs around one component of the binary rather than circumbi¬ 
nary discs. The specific parameters used here can be summarized 
as follows: The disc is initially planar and extends from R^^ - 0.1a 
to /?out = 0.35a with a surface density profile S = and 
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Figure 4. 3D surface rendering of the disc which was initially inclined at 60° to the binary plane with no warp. These snapshots are taken after 0, 2.5, 3.5 and 
5.5 binary orbits. The disc is viewed along the binary orbital plane and the arrow points in the direction of the companion in the snapshots. Four stages are 
shown: (i) the initial inclined disc, (ii) the disc is disrupted by the precession torque, (iii) the disc breaks into two distinct rings, but the outer ring is growing 
more eccentric and (iv) the outer ring and inner ring merge, causing enhanced dissipation and leaving behind an inclined eccentric disc. 


locally isothermal sound speed profile Cg = Cs,o{R/Rm)~^y where 
we have chosen /? = 3/2 and q = 3/4. This achieves a uniformly 
resolved disc with the shell-ave raged smoothing length per disc 
scale-height (h) IH ^ constant (iLodato & Pringle. 2007h . Eq and 
Cs,o are set by the disc mass, = 10“^M and the disc angular 
semi-thickness, H/R - 0.01 (at R - R\n) respectively. Initially the 
disc is composed of 4 million partieles, which for this setup gives 
(h) IH ^ 0.5. The simulations use a disc viscosity with Shakura & 
Sunyaev cr ^ 0.1 (wh ich requires artificial viscosity a ay = 1.9; cf. 
iLodato & Price 201^ and a quadratic artificial viscosity JSay = 2. 
We assume that the binary components, represented by two Newto¬ 
nian point masses with Mi - M 2 - 0.5M, aecrete any gas coming 
within a distance 0.05a of them, and so remove this gas from the 
computation. The simulations differ only in the initial inclination 
angle between the disc and the binary orbit. We perform our simu¬ 
lations for 6 - 10°, 30°, 45° and 60°. 

Fig. m shows the simulation with an initial inclination of 10°. 
Here the precession torque eaused by the eompanion is weak, so 


the disc evolves with a mild warp. We know from Eq. [6] that the 
strength of the precession torque is higher when 6 = 30°, and it has 
its maximum value when 6 - 45°. This agrees with Eq. [8] which 
shows that disc breaking is more likely when sin 2^ is high. We 
find strong disc breaking in our simulations with initial inclinations 
of 30° and 45°. Fig.|2| shows a simulation with an initial inclina¬ 
tion of 30°. Here the disc becomes signifieantly warped after a few 
orbits. Then the outer disc breaks off to form a distinct outer ring. 
Similarly, the disc with an initial inclination of 45° is disrupted by 
the strong precession torque and initially breaks into two distinet 
planes. Then a third ring is broken off, but quickly interacts with 
the outer ring. The two outermost rings merge after another ~10 
binary orbits, as shown in Eig. [3] The inner ring develops a strong 
warp, aligning somewhat towards the binary plane, while the outer 
ring remains highly inelined and preeessing. The 6 - 60° simu¬ 
lation evolves quite differently from those with smaller 6. This is 
shown in Eig. |4] The simulated disc appears as if it is about to tear 
after a few binary orbits, but the outer regions of the disc become 
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quite eccentric. The inner and outer disc interact strongly at peri- 
centre, causing enhanced dissipation (cf. Fig. [5]) and merging the 
two rings into a single eccentric disc. The remaining eccentric disc 
persists over the duration of the simulation, which is approximately 
50 binary orbits. 

The simplified criterion ([8]) derived in Section [2| predicts the 
breaking found in the simulations with inclination angle > 30°. 
However, it suggests that the disc with inclination angle 10° should 
also break with /^break ^ 0.23, within the disc outer radius (/?out = 

0.35). Instead of breaking the disc, we find that the 10° simulation 
aligns to the binary plane after ~ 20 binary orbits. This alignment ^ 
suggests that the vertical viscous torque ID) is dynamically impor¬ 
tant on short timescales. In this case, the simplifications made be¬ 
tween 0 and IS]) are not relevant. We also note that the outer disc 
radius in the 10° simula tion shrinks slightly to - 0.3a due to 
disc-binary resonances (lArtvmowicz & Lubow. 1994r) . 

The importance of the vertical viscous torque for the 10° sim¬ 
ulation can easily be shown by estimating the a 2 term in Eq. [7] 
Analysing this simulation we fi nd the warp am plitude grows to 
|j^| 0.1 which gives a 2 ~ 52 dOgilvie. 19991) . From these val¬ 

ues we find /^break ^ 0.41a by considering the vertical viscosity 
(see eq. [7]). The breaking radius predicted by the vertical viscous ^ 

torque exceeds the disc outer radius and this is in agreement with o 

the simulations. Therefore we can conclude that the simplified cri- | 

terion is not relevant for this case. However, for the 30° sim- g 
ulation we find that the warp amplitude grows to \4f\ 1.5 which g 

gives Qf 2 ~ 0.4 dOgilvie. 19991) . From these values we find that the ^ 

disc should break with /^break ^ 0.2a, within the disc outer radius, 
again in agreement with the simulation. Further we find that the 
30° simulation breaks at a minimum radius of 0.205a, whereas the 
45° simulation breaks at a minimum radius of 0.195a. From our 
estimate in Section [2] we expect /?break to differ between these two 
simulations by a factor of [sin 60/ sin90]"^/^ = 1.05, which is re¬ 
markably similar to the difference found in the simulations. These 
numbers give the smallest radius at which the disc was deemed to 
have broken, occurring at t = 19 in the 30° simulation and t = 16 in 
the 45° simulation. 

The accretion rate through tearing discs is generally signifi¬ 
cantly enhanced. The top panel of Fig. [5] shows the accretion rates 
for the simulations with 6 - 10°, 30° and 45°. The accretion rate 
is higher for the broken discs {6 = 30°,45°) than for the warped 
disc {6 - 10°). The disc with 6 = 60° produces highly variable 
accretion, as shown in the lower panel of Fig. O varying by ap¬ 
proximately three orders of magnitude. The high accretion rate for 
this simulation results from the enhanced dissipation between the 
inner disc and the eccentric outer regions. The remaining eccentric 
disc shows a nodding motion which produces the peaks in accretion 
rate seen in the bottom panel of Fig. [5] 

The nodding motion (oscillations in the disc tilt) are ac¬ 
companied by oscillations in the disc eccentricity. We attribute 
this to the Koza i-Lidov mechanism recently discovered by 
iMartin et al. f2014h to also act in fluid discs. The Kozai-Lidov 
mechanism induces antiphase oscillations in the orbital inclination 
and eccentricity of particles h ighly inclined to a binary compan¬ 
ion (iKozai. 196^1 5 dov. 19621) . iMartin et al. (20141) are the first to 
show that this process also occurs in a fluid disc. In our 60° sim¬ 
ulation the disc quickly becomes a narrow ring (through a strong 
interaction induced by tearing, which drives particles outside the 
ring into the accretion radius of the primary). This narrow ring 
then goes through strong Kozai-Lidov cycles. When the disc has 
its peak eccentricity the pericentre of the ring approaches the ac¬ 
cretion radius and creates the strong accretion rate shown in Fig. [5] 




time 


Figure 5. The accretion rates with time for different inclination angles of 
the disc. The accretion rate is calculated in time bins of half a binary orbit. 
The accretion rate is in code units of binary mass per binary dynamical time, 
and the time axis is in units of the binary dynamical time where a full orbit 
of the binary is 2n. 


IMartin et al. (2014h were the first to show that this behaviour oc¬ 
curs more generally for extended discs, and is therefore probably 
important in many astrophysical scenarios. 

The increase in accretion rate for the broken discs here is 
o nly a factor of a few , rather than the orders of magnitude found 
in iNixon et al. 120 13h . The main reas on for this is the w eaker pre¬ 
cession here (cf. eq.[T] compared with iNixon et al. 201^ . eq. 7). In 
the present case the inner disc evolves significantly faster, reach¬ 
ing a degree of alignment before tearing happens (see Fig.O. This 
makes the internal disc inclination angle smaller than 26. Another 
important difference is that resonances hold out a circumbinary 
disc, but in our case resonances do not slow accretion. We note 
that for smaller mas s ratios (Mp < M\ ) the resonances driving su¬ 
perhump behaviour (I Whitehurst. 1988h appear in the disc, and we 
will investigate this in a further paper. We note finally that signifi¬ 
cantly higher accretion rates might well appear in the disc for other 
parameters, e.g. if the viscosity driving alignment was smaller. 
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Figure 6. 3D surface rendering of a thick {HIR = 0.05) disc in the wave-like warp propagation regime {HIR > a = 0.01). The disc was initially inclined at 
45° to the binary plane and becomes slightly eccentric over the duration of the simulation. The left hand panel shows the initial conditions and then the middle 
and right hand panels show the disc after a quarter and a half of a precession period respectively, i.e. a precession of 90° and 180°. The disc is viewed along 
the binary orbital plane. 


3.1 Waves 


In this Section we report a single simulation performed with a 
thicker disc, H/R = 0.05 and lower dissipation, a = 0.01, 
but otherwise identical parameters to the 45° simulation above. 
For these parameters (h)/H = 0.17 and therefore a ay = 0.56. 
For this simulation, the differential precession induced by the bi¬ 
nary is expected to be communicated through pressure waves, 
propagating at a velocity ^ Cs/2, induced b y a warped disc 
ilPapaloizou & Lin. 1995Lubow & Qgilvie. 2000h . For the param¬ 
eters of the simulation, the wave travel time across the disc is 
tw ~ ^out/<^s(^out) = 1.3 binary orbits. However, the fastest pre¬ 
cession time induced in the disc is ^prec(^out) = 2;r/[Qp(/?out)] = 6.4 
binary orbits, and so we expect the disc to be able to communicate 
the precession efficiently throughout the disc by pressure waves. 
This leads to global precession of the disc as seen by previous in¬ 
vestigations (iLarwood et al.. 199^ . Fig.[6]shows half a precession 
period taking ^13 orbits, so the global precession period observed 
in the simulation is approximately 25 binary orbits. The predicted 
global precession time can be calculated by dividing the integral of 
the angular momentum by the integral of the torque across the disc. 
For the parameters chosen here, this gives 34 binary orbits, which is 
consistent with the observed period as the calculation is highly sen¬ 
sitive to the inner and outer disc radii and the surface density profile 
for our parameter choices and these can change somewhat during 
the simulation due to viscous spreading of the disc and accretion. 

The global precession observed in this simulation reiterates 
our comment above that the criterion derived in Section |2] should 
not be applied to scenarios for which the equations are invalid. 
This simulation is in the wavelike warp propagation regime with 
HjR > a. Therefore the internal disc communication is dominated 
by waves rather than viscosity. Therefore it is no surprise that this 
simulation does not agree with ([8l) as this was derived assuming 
that warps propagate diffusively. If we instead consider the internal 
disc communication due to waves, it is clear that the disc should 
not tear as the wave communication across the whole disc occurs 
faster than any local precession time in the disc. This is confirmed 
by our simulation shown in Fig.[6| 


4 DISCUSSION 

We have simulated misaligned accretion discs in a binary system 
to explore the process of disc tearing, where the precession torque 
induced in the disc can overwhelm its internal viscous communica¬ 
tion. Here a misaligned binary companion gravitationally drives the 
precession torque. We have shown that sufficiently thin and suffi¬ 
ciently inclined discs can break so that their outer rings precess 
effectively independently. This process can enhance the dissipation 
in the disc and promote stronger accretion on to the central object. 

In our analytical estimates in Section [2| we compared the pre¬ 
cession torque to the internal viscous torques arising from both az¬ 
imuthal and vertical shear. As it is not straightforward to calculate 
the vertical torque, which is strongly dependent on warp amplitude, 
we simplify the full criteria (Stoi). This simplification is relevant 
for moderate to large values of the disc viscosity parameter a and 
large inclination angles. We note that an exact calculation of a tear¬ 
ing criterion would require knowledge of dl/dR as a function of 
both position and time, so that the viscosity coefficients and hence 
the viscous torques can be calculated and compared to the preces¬ 
sion torque. For these reasons one has to perform three dimensional 
hydrodynamic simulations of this process to find out exactly what 
the disc does. The simplest approach to a criterion for disc tearing 
is given by i), which has the advantage of being readily calcu¬ 
lable, but is not applicable to all regions of parameter space. We 
shall explore this further with a focussed investigation in a future 
publication. 

We note that the criteria we have derived should not be used 
when the disc viscosity is smaller than the disc angular semithick- 
ness (a < H/R) as this allow s the efficient propagation of waves 
jPapaloizou & Pringle. 1983h and this distinct internal disc com¬ 
munication is not included in our analysis in Section |2] We have 
performed one simulation of such a disc with HIR = 0.05 and a = 
0.01 (Section [3TI ). If we were to naively apply the diffusive tearing 
criterion to this pressure dominated simulation we would ex¬ 
pect the disc to tear similarly to Figs [3] & |4] However, we instead 
find that the disc precesses as a solid body. This happens because 
the differential precession induced by the binary is communicated 
across the whole disc by pressure waves, which propagate at a ve - 
locity Uw ~ Cs/2 (jPapaloizou & Lin. 199 jPubow & Qgilvie. 2000l) . 
This leads to global precession of the disc, as seen in previous in- 
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vestigations llLarwood et al.. 19^ . We note that if we consider a 
criterion which tak es account of wave communication in the disc 
(iNixon et al.. 2013h . then we instead predict that the disc should 
not break, consistent with this simulation. 


5 CONCLUSIONS 

We have shown that tilted discs inside a binary are susceptible to 
tearing from the outside in, because of the gravitational torque from 
the companion star. If the disc inclination is small the disc warps 
(Fig. [T]), and if there is nothing to maintain the tilt, the disc even¬ 
tually aligns with the binary plane. For larger inclinations the disc 
can be torn, the outer ring breaking off and processing effectively 
independently (Figs.[2|&[3l). For some inclinations the eccentricity 
of the outer disc grows (Fig. 0. 

The behaviour we have discussed in this paper is relevant 
to a variety of astrophysical systems, for example X-ray bi¬ 
naries, w here the disc plane ni a y be tilted by radiatio n warp¬ 
ing (e.g. IWiiers & Pringle 19^ : fOgilvie & Dubus 200lh . super- 
massive black hole binaries, where accretion of misaligned gas 
can create effectively random inclinations (e.g. iNixon et al.. 2013h 
and protostellar binaries, where a disc may be misaligned by 
a variety of effects such as binary captur e/exchange, accre- 
tion after binary forma tion (iBate et al.. 201C)I) and stellar flybys 
(INixon & Pringle. 20101) . 
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